Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 10 August 2012 



(MN KTeX style file v2.2) 



Influence of XUV radiation on P v ionization fraction in hot star 
winds 



. Jin Krticka and Jin Kubat 

■ ^Dstav teoreticke fyziky a astrofyziky, Masarykova univerzita, CZ-611 37 Brno, Czech Republic 
' ^Astronomicky ustav AV CR, Fricova 298, CZ-251 65 Ondfejov, Czech Republic 



Received 



ABSTRACT 

Different diagnostics of hot star wind mass-loss rates provide results that are difficult to 
reconcile with each other. The widely accepted presence of clumping in hot star winds implies 
a significant reduction of observational mass-loss rate estimates from diagnostics that depend 
on the square of the density. Moreover, the ultraviolet P V resonance lines indicate a possible 
need for even stronger reduction of hot star mass-loss rates, provided that P V is a dominant 
ionization stage of phosphorus at least in some hot stars. The latter assumption is challenged 
by a possible presence of the XUV radiation. 

Here we study the influence of the XUV radiation on the P V ionization fraction in the 
hot star winds. By a detailed solution of the hydrodynamical, radiative transfer, and statistical 
equilibrium equations we confirm that sufficiently strong XUV radiation source may decrease 
the P V ionization fraction, possibly depreciating the P V lines as a reliable mass-loss rate indi- 
cator. On the other hand, the XUV radiation influences also the ionization fraction of heavier 
ions that drive the wind, leading to a decrease of the wind terminal velocity. Consequently, 
we conclude that the XUV radiation alone can not bring theory and observations in accord. 

We fit our predicted wind mass-loss rates by a suitable formula and compare the results 
with the observational mass-loss rate diagnostics. We show that for supergiants and giants the 
theoretical predictions do not contradict the mass-loss rate estimates based on X-ray line pro- 
files or density squared diagnostics. On the other hand, for main-sequence stars the predicted 
mass-loss rates are still significantly higher than that inferred from P V or X-ray lines. This 
indicates that the "weak wind problem" recently detected in low-luminosity main-sequence 
stars may occur to some extent also for the stars with higher luminosity. 

Key words: stars: winds, outflows - stars: mass-loss - stars: early-type - hydrodynamics - 
X-rays: stars 



1 INTRODUCTION 

Mass-loss plays an important role in the massive star evolution. 
Most of the mass of massive stars is lost during their evolution from 
zero-age main sequence to the final remnant. Different processes 
contribute to the mass loss in individual evolutionary phases. These 
processes include line-driven winds during hot evolution ary stages 
jPuls et al.ll2008ah. decretion disks in fast-rotating stars I ILee et al.l 



iQkazak illiqo'lh. LBV-type of explosions in hot supergiants 



dSmith & Owockil 1200^ . dust-driven winds i n cool supergiants 
(Woitke 2006), and a final supernova explosion iUmeda & Nomotol 
[2OO8). Consequently, estimates of amount of mass lost per unit of 
time (mass-loss rate) as a functions of stellar parameters belong to 
one of the most important ingredients of evolutionary models. 

Unfortunately, the uncertainties in modern mass-loss rate de- 
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terminations significantly affect the evolutionary models of mas- 
sive stars. In the case of the line-driven wind of hot stars, these 
uncertainties seem to be most ly connected with th e occurrence of 
small-scale inhomogeneities jHamann et al.l I2OO8I) . The inhomo- 
geneities are typically divided into three differents groups (micro- 
clumping, porosity, and vorosity) according to their influence on 
the spectral features, although they may be caused by the same 
structure observed in different wavelengths. Microclumping (also 
frequently referred to as clumping), that accounts for the enhanced 
density in the optically thin inhomog eneities, can be most easily 
incorporated in t he wind models (e.g ., iHamann & Grafene3l2004l : 



IPuIs at all [200^ : iKrticka et alj 12003) . Microclumping affects the 
ionization equilibrium via enhanced recombina tion, consequently it 
influences the radiative transfer only indirectly dAbbott et al.ll98lh . 
On the other hand, porosity (also referred to as macroclumping) 
accounts for the nonnegligible optical depth of inhomogeneities 
(which may become optically thick), and directly infuences the 
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radiative transfer jOskinova et alj |2007| ; ISundqvist et alj llOllI : 
ISurlanetal.ll2012ah . Vorosity affects the line profiles, and it is 
connecte d with the diffe rent Doppler shifts of individual inhomo- 
seneities i lOwockill2008h . From the point of view of mass-loss rate 
predictions, the inhom ogeneities may affect the ionization fract ions 
of wind driving ions jKrticka et al.ll2008l ; iMuiires et alj|201 ih , or 
they may lead to the decrease of the wind mass-loss rate due to the 
base turbulence ( lLucvll2007l : lKrticka & Kubi3l2O10l) . 

There could be a simple observational solution of the mass- 
loss rate determination problem: to find such observational charac- 
teristic that is not affected by the wind inhomogeneities. There are 
two potential candidates for such convenient observables: the X- 
ray radiation (either line profiles or the continuum flux distribution, 
MacFarlaneetal."l99l"; 'Owo cki & CohenI 1200 ll : llgnace & GavlevI 
2002; Cohen et al. 20 11 a), and unsaturated resonance line profiles. 
The latter case is fulfilled for trace elements , and especially in- 
teresting is the Pv i on tCrowther et alj|2002l : iMassa et al.ll2003l : 
lFullertonetn]|2006h . However, even these characteristics face 
some problems. The opacity in the X-ray domain scales mostly 
linearly with the density, consequently its effect on the X-ray diag- 
nostics should be in principle possible to model in a straightforward 
way. Indeed, the observations of mostly symmetric X-ray line pro- 
files, that are not strongly affected by absorption, indicate low wind 
mass-loss rates jWaldron & Cassinellil 1200 ll ; ICohen et alj |2010| ; 
iGagne et al.l2oT l|). These results are, however, chal lenged by a pos- 
sible effect of the porosity on the X-ray opacity jFeldmeier et al.l 
l2003l ; lOskinova et alj|200q). Note that there is not a general con - 
sensus on this problem iOskinova et alj2012l ; IOwocki et al.ll2012h . 

Here we concentrate mostly on the other promising observa- 
tional characteristic, i.e., on the Pv resonance lines. The weak- 
ness of observ ed Pv lines also ind icates low mass-loss rates of 
hot star winds jFullerton et al.l2006l) . However, the ioniz ation frac- 
tion of P V may be modifie d by the effect of clumping ( Puis et all 
l2008bl ; iKrticka et al.l I2OO8I) . Additional changes of the Pv ion- 
ization fraction may be caused by the influence of X-rays. Al- 
though our previous calculations indicate d that P V ionization frac - 
tion is not stro ngly affected by X-rays jKrticka & Kubatll2009h . 
IWaldron & Ca ssinelli (2010) argued that the extreme ultraviolet ra- 
diation (hereafter XUV) may affect P V io nization fractions. Since 
the work of iWaldron & Cassinellil Bm(X) was based on rather 
simplified ioniza t ion es timates and the calculations presented in 
iKrticka & Ku"bi3 ( l2009h were done without enhanced XUV radi- 
ation, we decided to fill this gap and we apply our NLTE wind 
models to study the influence of the XUV radiation (parameterized 
in a convenient way) on the P V ionization fraction. The XUV re- 
gion is defined here as the energy from interval 54.4 eV (He II edge) 
to 124eV. 



2 WIND MODELS 

For o ur calculations we use NUTE wind models of lKrticka & KubiJ 
with a comoving frame (CMF) line force. Our models as- 
sume stationary, and spherically symmetric wind flow. They enable 
us to selfconsistently predict wind structure just from the stellar 
parameters (the effective temperature, mass, radius, and chemical 
composition). The line radiative force is calculated directly by sum- 
ming the contribution from individual atomic transitions, i.e. we do 
not use the CAK line force parameters. 

The ionization and excitation state of the considered elements 
is derived from the statistical equilibrium (NLTE) equations. Ionic 
models are either adopted from the TLUSTY grid of model stellar 



Table 1. Stellai' parameters of the model grid 



model Toft 
[K] 



[Rq] 



M 

[Mo] 



M 

[M0 year^-"-] 



mam 
sequence 



300-5 
325-5 
350-5 
375-5 
400-5 
425-5 



30000 
32 500 
35 000 
37 500 
40000 
42 500 



6.6 
7.4 
8.3 
9.4 
10.7 
12.2 



12.9 
16.4 
20.9 
26.8 
34.6 
45.0 



9.5 X 10-9 
9.8 X 10-9 

4.6 X 10-* 
1.8 X 10-^ 
6.4 X 10-^ 
1.2 X 10-*^ 



giants 



300-3 
325-3 
350-3 
375-3 
400-3 
425-3 



30000 
32 500 
35 000 
37 500 
40000 
42 500 



13.1 
134 
13.9 
14.4 
15.0 
15.6 



19.3 
22.8 
27.2 
32.5 
39.2 
47.4 



8.7 X 10-* 

1.8 X 10-^ 
4.5 X 10-^ 
1.0 X 10-s 

1.8 X 10-'^ 

2.9 X 10-*^ 



supergiants 



300-1 
325-1 
350-1 
375-1 
400-1 
425-1 



30000 
32 500 
35 000 
37 500 
40000 
42 500 



22.4 
214 
20.5 
19.8 
19.1 
18.5 



28.8 
34.0 
40.4 
48.3 
58.1 
70.3 



4.7 X 10-^ 
7.7 X 10-^ 
1.3 X 10-^ 
2.3 X 10-'^ 
3.0 X 10-'^ 
3.7 X 10-^ 



atmospheres dUanz & Hubenvll2003Ll2007h or a re created by us us- 
ing the data fro m the Opacity and Iron Projects /Seaton et a l.ll99j ; 
[ifummer et al. 1993h. For p hosphorus we employed data described 
by Pauldrach et al. (l200lh . Auger photoioniz ation cross sections 
from i ndividual inner- shells were taken from IVemer & Yakovlevi 
( I1995I see alsolVemeretal.l | l993h . and Auger yields were taken 
from lKaastra & Mewd ( 119931)" The emergent surface flux is taken 
fr om H-He spher ically symmetric NLTE model stellar atmospheres 
of iKubaJ ( I2OOI and references therein). For our wind calcula- 
tions w e assume a solar chemical composition after lAsplund et all 
( I2OO9I) . 

The radiative force is calculated using the solution of 
the spherically sym metric CMF radiative transfer equation 
dMihalas et"^ll975h . The corresponding line data were extracted 
in 2002 from the VALD database (Piskunov et al. ll995L Kupka et 
al. ll999l) . The radiative cooling and heating terms are deriv ed us- 
ing the electron thermal balance method (Kubat et al., 1999). For 
the calculation of the radiative force and the radiative cooling and 
heating terms we use occupation numbers derived from the statisti- 
cal equilibrium equations. The hydrodynamical equations, i.e., the 
continuity equation, equation of motion with the CMF line force, 
and the energy equation with radiative heating and cooling included 
are solved iteratively to obtain the wind density, velocity, and tem- 
perature structu re. The wind mass-loss rate is derived from the crit- 
ical condition dCastor. Abbott & Klein|[l973 ) generalized for the 
case of CMF line force. The derived mass-loss rate corresponds 
to the m aximum one for which smooth transonic solution can be 
obtained dPoe et al.lll99lh . 

For our study we selected O star model grid with the effective 
temperatures in the range 30 000 — 42 500 K. The parameters for 
the stars with give n effective temperatur es were obtained using re- 
lations derived by iMartins et alj d2005a[) for main-sequence stars, 
giants, and supergiants (see Table[T|l. 

Our new mode ls predict slightly lower mass-loss rate than 
our older models ( Krticka et al.l l2009 h due to inclusion of line 
overlaps via the solutio n of the CMF radiative transfer equation 
dKrticka & Kuball I2OIOI1 . The mass-loss rate predictions for all 
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models are also listed in the Table[T] We fitted these mass-loss rate 
predictions as 



log 



M 



W-^ Mq year 



= {a + ail)\og ^ 



L 



106 Lr. 



+ b + bil, 
(1) 

where I is the luminosity class (i.e., 1 for supergiants, 3 for giants 
and 5 for main-sequence stars), and 



a = 2.040, 
b = 0.552, 



ai = -0.018, 
bi = 0.068. 



(2) 



We compared for mula Eq. (TJ with mass-loss rate predictions of 
IVink et al. I ( I2OOII) calculated for the mass-fraction of heavier ele- 
ments Z = 0.0134 (note that the sola r mass-fraction of heavier 
elements assumed bv|Vinket alj ( 1200 ih is different, Z = 0.0194, 
I Anders & Grevessdl 19891) 7 On average our models predict slightly 
lower mass-loss rates by a factor of about 1.7. Formula Eq. (TJ 
shows an excellent agreement with lPauldrach et all ( 1201 2h models 
A"'"/A~ of Pup, while the predictions of our formula are by a fac- 
tor of about 1.6 lower than the mass-loss rate of models D^/D~. 



3 ADDITIONAL X-RAY/XUV RADIATION SOURCE 

We include an additional source of X-ray/XUV radiation into our 
wind models. For this purpose we use the X-ray emissivity -qx (r, u) 
dKrti cka et al. 2009, Eq. 11), which was derived using the n umer- 
ical simulations of wind instability dFeldmeier et alj[l997l) . This 
additional high energy emission starts at the radius 2 /?* , and in- 
tegrates the emission from the gas with different shock tempera- 
tures as derived from hydrodynamical simulations. In our approach 
the X-ray emission lines, which mostly contribute to the high en- 
ergy emission in hot stars, are not treated individually, but they are 
summed over selected wavelengths. This approach is fully satisfac- 
tory because most of these lines are optically thin in the cold wind, 
consequently the X-ray emission lines influence the wind ioniza- 
tion equilibrium via A uger and direct photoionization only (e.g., 
lMacFarlaneetalJl99lh . 

The original expression for the X-ray emissivity rjx (r, u) is 
modified to test the influence of the XUV radiation on the wind 
ionization. We introduce a nondimensional free parameter /3xuv 
that scales the original X-ray emissivity rjxij-, v) in the XUV re- 
gion. This region includes frequencies /^xuv < < i^xu\^ where 
VxiN is the Hell ionization threshold (54.4 eV), and we selected 
fxuv = 3 X 10^^ s~\ which corresponds to 124eV. The modified 
X-ray emissivity included in our models is therefore given by 

nx{r, v) = <^ /3xuvr?x(r, v), z.^]? <v< J'xuv , 0) 
[ r?x(r,i/), v>vxv^- 

For /3xuv = the additional XUV emission is not considered at all 
and only additional X-ray ionization source is included, whereas 
/3xuv > 1 corresponds to enhanced source of the XUV emission. 



4 C V IONIZATION FRACTION AND ADDITIONAL 
IONIZATION SOURCES 

The presence of lines of ions with higher degree of ionization ( for 
example C V) is often used (e.g., bv lWaldron & Cassinellil2()lol) as 
an argument for the existence of additional XUV or X-ray ioniza- 
tion source. Although it is generally true that additional ionization 
source shifts the degree of ionization, due to the complexity of the 



processes involved we can not claim that the presence of a particu- 
lar ion is caused exclusively by radiation at a chosen frequency. 

As an example let us study the ionization ratio of CiV and 
C V. Assuming that the population of the ground levels of these ions 
dominate, i.e. that we can neglect population of excited states, the 
ionization balance between the C IV ion (A'^4 is its number density) 
and C V ion (number density N^) fo llows from the equations of 
statistical equilibrium dMihalasll 19781) as 



iV4i?45 - iV5i?54 = 0, 



(4) 



where we took into account only the radiative ionization and re- 
combination (collisional transitions are neglected). The radiative 
ionization rate 



i?45 = 47r 



hv 



J{u) du, 



and the radiative recombination rate 



R54 



■ 4tv 



hv 



2hv^ 



J{u) 



(5a) 



e'T^ Av, (5b) 



where asterisk denotes LTE values and 04^1, is the photoionization 
cross-section. Replacing the integrals in Eq. (jsj with values at the 
ionization frequency 1^4 and taking into account that for a consid- 
ered spectral range 2hv''' jc? ^ J{y\ we derive from Eq. (|4j 



e kT 



(6) 



Using the Saha-Boltzmann equation, the fraction ( / A'4 ) * can be 
eliminated and the latter equation can be further simplified (assum- 
ing unity ionic partition function) to 

,(2™efcT)^/^:^, (7) 



N4 



where A'c is the electron density. From this equation it seems that 
indeed the ionization ratio is directly proportional to the mean ra- 
diation intensity J at a given ionization frequency 1/4. Using values 
appropriate for the model 350-1 at roughly 2.3 i?,, we obtain in 
absence of additional ionization sources for Vi = 1.6 X 10^® Hz, 
r = 22 000 K, iVe = 3 X 10^cm~^ and J{vi) = 8 x 
10^^^ erg cm^^ s^^ Hz from Eq. (|7]l the ionization ratio N5/N4 ~ 
7 x 10~*. This indicates low ionization fraction of C V there. How- 
ever, from our full NLTE models, which consider reliable model 
ions, we obtain for the same location N5 /N4, ~ 0.2. 

The reason is that statistical equilibrium equations are quite 
complex in hot star winds and their oversimplification using handy 
equations like Eq. Q may lead to incorrect results. In the particular 
case of C V ionization fraction, the most important ionization pro- 
cess is not that from the ground level, but from the less populated 
upper levels, which are closely coupled with the ground level by a 
strong bound-bound transitions. This information can be obtained 
only solving the equations of statistical equilibrium. 

Consequently, a care has to be taken when making the conclu- 
sions about the existence of additional ionization source just from 
the observations of C V li nes. The same comment is valid also for 
Nv lines ( Pauldrachlll987i ; lKrticka & KubafcoOSi) . 



5 INFLUENCE OF XUV RADIATION ON P V 
IONIZATION FRACTION AND WIND MODELS 

Let us now turn our attention to Pv. Without additional XUV 
sources, for /?xuv ~ our models agree with the conclusions of 
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Figure 1. Ionization fraction of individual pliosphorus ions in the model 
400-1 with different XUV sources 
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Figure 2. Comparison of predicted product of averaged P V ionization frac- 
tion and the wind mass-loss rate g(Pv)M calculated with /3xuv = 1 
(empty circle s) and feuv = \/T0 (full circles) with results derived from 
observations jpuUerton et al.l2006l errorbars and upper limits) as a function 
of stellar effective temperature. The results for different luminosity classes 
are separated into individual graphs. 



lKrticka&Kuball l l2009 '). according to which the X-rays do not sig- 
nificantly influence the ionization fraction of P V. Only the ioniza- 
tion fraction of P VII (and P VI in cooler stars) is significantly in- 
fluenced by the X-rays, while P V stays to be one of the dominant 
ionization stages. 

The inclusion of XUV radiation (with /3xuv = 1, see Fig. \^ 
leads to a slight reduction of the P V ionization fraction to roughly 
60% of its value calculated without XUV radiation (/3xuv = 0). 
This indicates that a stronger XUV source may lead to a signifi- 
cant reduction of P V ionization fraction. Indeed, the P V ionization 
fraction is significantly reduced in the models with /3xuv ~ VT^, 
as shown in Fig. |2] Here we plot the product of P V ionization 
fraction and wind mass-loss rate g(Pv)Af averaged over radii 
2 < r/J?« < min(5,a;max) with a;max being the radius of the 
outer model boundary (in units of _R«). 

The comparison of predicte d Pv ionization frac tions and 
those derived from observat ions by Fullerton et al.l ( 
supports the suggestion of IWaldron & Cassinellil I 



2006) in Fig. I2l 



20101) that en- 



hanced source of XUV radiation leads to a reduction of the P V 
ionization fraction. There is a good agreement between the product 
g(P v)Af for the para meter /3xuv = vTO and the observational one 
(determined by Fuller ton et alj20()^ for supergiants and giants. On 
the other hand, the results for main-sequence stars are ambiguous. 
While for main-sequence stars with detected P v line there is a good 
agreement with observations even without additional X-ray/XUV 
source, for stars with only upper limit of q(P v)M available even a 
very strong X-ray/XUV source does not bring the predictions and 
observations into agreement. Suc h disagreement can be a signature 
of a "w eak wind problem" (e.g. iBouret et alj|2003l ; iMartins et al.l 
l2005d) . 

IWaldron & Cassinelll ( 1201 Oh argued that due to a unique dis- 
tribution of XUV radiation, where the strongest XUV emission 
lines have energies lower than the S V edge (72.7 eV), the XUV 
radiation affects the ionization fraction of P V and not that of S V. 
This effect might not be fully included in our models, where the 
XUV emission lines are summed over a corresponding wavelength 
region. This ensures that the total line emissivity is properly taken 
into account, but some subtle effects of line distribution may be 
missing. To test the effect of XUV line distribution, we calculated 
additional models with XUV emission present only for energies 
lower than 72 eV. These models confirm that with a convenient 
distribution of XUV emission the influence of XUV radiation on 
P v ionization fraction is more significant than on S V (roughly by 
a factor of 2). However, even the ionization fraction of S V is af- 
fected by XUV radiation with energies lower than the ionization 
energy of S V due to the ionization from higher levels of S V, espe- 
cially from the relatively strongly populated second level 3s3p "^P 
with a ionization energy 62.4 eV. This is a similar situation to C V 
ionization fraction discussed in Sect.|4] 

However, despite the promising results derived for P V in gi- 
ants and supergiants, a detailed inspection of our models (which 
consistently include all possible driving ions) shows that there ex- 
ists one additional (and natural) effect of the XUV radiation. Phos- 
phorus is a trace element and changes of its ionization balance 
have only negligible effects on the radiation force. However, also 
many ionization states of non-trace elements, which are the wind 
drivers, are depopulated by XUV radiation. These ions include 
Civ (ionization energy £;,o„ = 64.5 eV), Niv (S,o„ = 77.5 eV), 
O III (Sion = 54.9 eV), O IV (Ei„„ = 77 A eV), and Si IV (£,„„ = 
45.1 eV). This causes drastic changes in the line force accelerating 
the whole wind. Consequently, the wind becomes overionized and 
the radiation driving by higher ionization states becomes inefficient 
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Figure 3. Influence of XUV radiation on the wind radial velocity in tlie 
model 400-1. 



iLeutenegger et alJllOOd : ICohen et al.ll201 la l). By contrast, the on- 
set of XUV emission closer to the star suggested by observations 
leads to even lower wind terminal velocity, which does not agree 
with observations. 



6 DISCUSSION: WIND MASS-LOSS RATES - THEORY 
MEETS OBSERVATIONS? 

The question of correct mass-loss rate determination is the most 
important one of any wind theory and observation. Despite a sig- 
nificant progress in line driven wind theory and especially in our 
under standing of the radiative transfer in structu red winds in recent 
years l lSundgvist et al.l201 iMSurlan et al.l2012al) . the answer to this 
question is not unambiguously solved. Concerning the P V ioniza- 
tion fractions studied here, a reduction of mass-loss rate predictions 
by a factor of about 3 is necessary to bring the observation of giants 
and supergiants and theory into agreement (assuming /3xuv = 1 not 
to violate the observed terminal velocities). 
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Figure 4. Comparison of predicted ratio of the wind terminal velocity and 
escape velocity calculated with /3xuv = 1 (empty circles) and /3xuv = 
\/lO (full circles) with results derived from observed Vaa jLamers et alj 
[199^ . 



due to their insufficient line opacity. This causes significant lower- 
ing of the radiative force. Consequently, the wind overionization 
leads to wind stagnation that starts roughly at the radius where the 
XUV radiation is switched on in our models, i.e., at r = 2R, (see 
Fig.O. In our calculations models of giants and supergiants of all 
spectral types and main sequence stars with Tcff > 35 000 K are 
subject to this effect. 

Such stagnation would be observationally manifested by a de- 
crease of the wind terminal velocity. In Fig. |4] we plot the ratio 
of the wind terminal velocity to the escape velocity Woo/vesc from 
our models in comparison with observational results. For a weak 
additional XUV source Pxuv = 1 the predicted Vpo /ve sc nicely re- 
produces the observed results jKrticka & Kubi3 l200^. However, 
for a stronger XUV source /3xuv = VW that reproduces the ob- 
served Pv ionization fractions, the predicted Voo/vex is signifi- 
cantly lower than that based on observational results. This disagree- 
ment could be in principle partially compensated in models with 
an onset of additional XUV emission farther from the s tellar sur- 
face, however this is supported neithe r by models iFeldmeier et al.l 
Il997h nor by the observations (e.g., IWaldron & Cassinellill200ll : 



6.1 P V ionization fractions with XUV emission and clumping 

IPuIs et ai] ( l2008bl) proposed that the probl em of weak Pv lines 
is caused by mic roclumping (see also ICrowther et alj |2002| ; 
iKrticka et al High density inside the clumps favors the re- 

combination leading to decrease of Pv ionization fraction. We 
tested if the combination of microclumping and additional XUV 
source does not bring observations and theory into agreement. For 
this purpose we calculated other models, in which we both allowed 
for additional XUV emission and took into account the influence 
of higher wind den sity in the clumps on the NLTE equations (as in 
iKrticka et alj2008h . Within this microclumping approach the inho- 
mogeneities directly affect the ionization equilibrium only and do 
not influence the radiative transfer due to porosity. We assume that 
the clumping starts above the critical point at the same radius as 
the additional X-ray emission, i.e., at 2 J?,. In this ca se the clump- 
ing does not affect the predicted mass-loss rates (cf. IKrticka et al.l 
2008). 

Our models show that clumping does not bring observations 
and theory into agreement. Clumping has an opposite effect than 
additional XUV source decreasing the wind ionization. Therefore, 
with microclumping (and XUV emission) the phosphorus ion- 
ization fractions disagree with values derived from observations, 
while the higher radiative force (due to lower ionization) provides 
better agreement with observed terminal velocities. We were un- 
able to find such combination of parameters (describing additional 
XUV emission and microclumping) that would provide both phos- 
phorus ionization fractions and terminal velocities in agreement 
with observations. On the other hand, the porosity, i.e., the ef- 
fect of wind structure on the line formatio n jOskinova et al.ll2007l ; 
ISundgvist et"ai]|201 ll; Isurlan et al.ll2012bl) is the more promising 
one for the explanation of the remaining discrepancy between the- 
oretical and observed P V line profiles. 

6.2 Mass-loss rates from X-ray diagnostics 

The situation is a bit different in the case of the mass-loss rate 
derived from X-ray line profiles than in the case of P V line pro- 
files The shape of X -ray line profiles may also imply low mass- 
loss rates in hot stars jWaldron & Cassinelli ^2001). Contrary to P V 
line profiles, the microclumping does not affect the shape of X- 
ray line profiles. However, the influence of porosity is still debated 
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Table 2. Wind mass-loss rates estimated for individual stars from X-ray diagnostics compared with our prediction after Eq. (T). 



Star 


Sp. type 


log(L/lL0) 


A/[lO"^M0year-i] 
X-ray diagnostics prediction 


Source (luminosities and X-ray 
mass-loss rates) 


HD93129A 


02If 


6.17 


6.8 ± 2.5 


9.2 


Cohen et al. (2011a), Repolust et al. (2004) 


CPup 


04If 


5.86 


3.5 ±0.3 


2.2 


Cohen et al. (2010), Puis at al. (2006) 


HD 93250 


04IIIf 


5.95 


1.4 ± 0.5 


6.2 


Gaeneetal (2011) 


9Sgr 


04V 


5.67 


0.34 


1.8 


Cohen et al. (201 lb), Martins et al. (2005a) 



dOskinova et"ai]|2012l ; IOwocki et al.ll2012t) . Moreover, these deter- 
minations involve some simplifying approximations, for example, 
the constant ionization structure of the ambient cool wind. 

In Table [2] we compare the mass-loss rates calculated us- 
ing formula Eq. lUll with those derived from X-ray line profiles 
and from continuum X-ray absorption in the case of HD 93250. 
The comparison may be biased, because HD 93250 and 9 Sgr are 
binaries and show peculiar features, including enhanced o r non- 
thermal radio emission or colidirig winds ([Le itherer et aP 1 19951: 
Rauw et al.ll2002l : ISana et aljUol ll : lOagne et ^.11201 II : iRauw et al.l 
2OI2I) . Note also that the mass-loss rate in the case of C, Pup was 
deriv ed from observatio ns assuming non-solar chemical composi- 
tion iCohen et allboiOl) , while our models assume solar chemical 
composition. 

The results shown in Table |2] show that there is a good agree- 
ment between wind mass-loss rates derived from observations 
and theory for the two supergiants, while for the other two non- 
supergiant stars predicted rates overestimate the observed ones by 
a factor of about 5. This result may reflect a similar dichotomy be- 
tween supergiant and main-sequence mass-loss rates found from 
P V ionization fractions (see Fig.|2]l. 



CO 
0) 



CO 
0) 



10-' 



10'- 



10-^ 



10' 



10'- 



10'' 



10' 



— 1 ' 1 ' 1 ' r- 



+ Puis et al. (2006) 
prediction 



luminosity class I 

H ' 1 ' 1 ' 1 ' 1 ' — 



+ Puis et al. (2006) 
prediction 



luminosity class II 



1.2 5.4 5.6 5.8 6.0 6.2 6.4 

log (/VLg) 



6.3 p diagnostics as upper limit to the mass-loss rates 

The Ha emission line and infrared (or radio) continua belong to 
traditional mass-loss rate indicators. However, the amount of Ha 
emission or infrared excess is not di rectly proportional to the wind 
density p, but to its square (? (e.g.. |Puls et alj|2008lj) . The struc- 
tured (clumped) wind with lower mass-loss rate may mimic spec- 
tral features of the homogeneous wind with higher mass-loss rate. 
Consequently, the Ha emission line or infr ared continua ma y pro- 
vide only upper limits to the mass-loss rate jPuls at al.l2()0^ , since 
they give values which correspond to a homogeneous wind. The 
predicted mass-loss rates should be lower than that derived from 
(? diagnostics (i.e.. Ha emission line or infrared continua). 

In Fig.|5]we plot the predicted mass-loss rates calculated using 
Eq. lO as a function of stellar luminosity in compariso n with ob- 
servat ional upper limits derived from (? diagnostics bv' Puls at al.l 
(l2006h . The luminosity of stars from this observa tional sarnple wa s 
derived using the stellar parameters taken from lPuls at al.l j2006l) . 
For supergiants with higher luminosity log(L/LQ) > 5.8 and for 
giants the mass-loss rates derived from observations are on aver- 
age only slightly higher than predicted ones. This indicates only 
medium level of clumping in these stars. On the other hand, the 
mass-loss rates derived from observations are significantly higher 
than predictions for low-luminosity supergiants \o%(L I Lq) < 5.8, 
pointing to a more significant level of clumping in their winds. 
However, in general the predicted mass-loss rates are not higher 
than those derived from the f? diagnostics. Consequently, even in 
this case the predictions do not contradict the observations. 



Figure 5. Comparison o f upper limits to the mass-loss rates from (?" diag- 
nostics jPuls at alJbood) with our predicted dependence of mass-loss rates 
Eq. Q as a function of stellar luminosity. 



7 CONCLUSIONS 

We tested the influence of XUV radiation on P V ionization frac- 
tions in hot star winds. Using hydro dynamical NUTE wind mod- 
elling we confirmed the conclusion of lWaldron & CassinellildioiOl) 
that the XUV radiation may decrease the Pv ionization fraction. 
As the result of various ionization and recombination processes, 
the P V ionization fraction in hot star winds does not come close to 
the value of 1, implying that the role of the Pv lines as mass-loss 
rate indicators in hot star winds is not so straightforward. On the 
other hand, the large amount of XUV radiation necessary to signif- 
icantly lower the Pv ionization fraction to the values implied by 
observations leads to a decrease of the wind terminal velocity due 
to inefficient line driving. This contradicts the observations. More- 
over, the wind clumping has an opposite effect than the XUV emis- 
sion decreasing the wind ionization. Consequently, it is unlikely 
that the effect of lowering the P V ionization fraction by XUV radi- 
ation alone can bring the theory and observations in accord. 

We also provide a useful mass-loss rate formula and compare 
it with other mass-loss rate diagnostics. We show that for super- 
giants and giants this formula passes two important observational 
tests against mass-loss rates derived from X-ray line profiles and 
(? diagnostics. This supports the reliability of mass-loss rates pre- 
dictions derived from modern wind codes for luminous hot stars. 

For main-sequence stars the predicted mass-loss rates are sig- 
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nificantly higher than those inferred from the Pv or X-ray hnes. 
This may indicate that the "weak wind problem" recently detected 
in low-luminosity main-sequence stars occurs to some extent also 
for the stars with higher luminosity. While the explanation of this 
"weak wind problem" for low-l uminosity main-sequence stars is at 
hand ( t oo large cooling length, Lucy & White 1980. Martins e t al] 
l2005bi , ICohenetalJ boost iKrticka & Kub at 2009. Lucy 201 3 . or 
influence of X-rays on the mass-loss rate, lDrew et al.il 19941) such 
explanation valid for all main-sequence stars with any luminosity 
is currently missing. 
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